Peculiar Nature of Snake States in Graphene 
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We study the dynamics of the electrons in a non-uniform magnetic field applied perpendicular to 
a graphene sheet in the low energy limit when the excitation states can be described by a Dirac type 
Hamiltonian. We show that as compared to the two-dimensional electron gas (2DEG) snake states 
in graphene exibit peculiar properties related to the underlying dynamics of the Dirac fermions. 
The current carried by snake states is locally uncompensated even if the Fermi energy lies between 
the first non-zero energy Landau levels of the conduction and valence bands. The nature of these 
states is studied by calculating the current density distribution. It is shown that besides the snake 
states in finite samples surface states also exist. 
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PACS numbers: 81.05.Uw, 73.21.-b, 73.63.Nm, 73.43.Cd 



In recent experiments on graphene new transport phe- 
nomena resulting from the massless Dirac fermion type 
excitations have been observed generatinghuge interest 
both experimentally and theoretically [H,ll, 01] . For re- 
views on graphene see Refs. 0, H, H and a special issue in 
Ref. 0. 

While for 2DEG a number of experimental and theo- 
retical works have been devoted to study the excitation 
spectrum of electrons and their transport properties in 
non-uniform magnetic fields, only a little is known about 
graphene in this case. For example, in 2DEG a special 
state called snake state exists in non-uniform magnetic 
field at the boundary where the direction of the mag- 
netic field changes. The snake states in 2DEG were first 
studied theoretically by Miiller [8| and it inspired numer- 
ous theoretical and experimental works (see, eg, Ref. [§| 
and references therein). The effect of an inhomogeneous 
magnetic field on electrons has been much less investi- 
gated in graphene. Martino et al. have demonstrated 
that the massless Dirac electron can be confined in inho- 
mogeneous magnetic field jlOj. Tahir and Sabeeh have 
studied the magneto-conductivity of graphene in a spa- 
tially modulated magnetic field and shown that the am- 
plitudes of the Weiss oscillation are larger than that in 
2DEG 11]. The low energy electronic bands have been 
studied by Guinea et al. [l2| taking into account the non- 
homogeneous effective gauge field due to the ripples of the 
graphene sheet. However, snake states have been stud- 
ied only in carbon nanotubes very recently by Nemec and 
Cuniberti [13j . 

In this work we study the electronic properties of 
graphene in a non-uniform magnetic field as shown in 
Fig. Q] We show the existance of snake states exibiting 
peculiar behaviour at low energy. In particular, we find 
that these snake states are localized in the zero magnetic 
field region and carry current which is compensated at 
the edges of the sample far from the central region. 

To study the nature of the snake states in graphene 
we consider a simple step like, non-uniform magnetic 



field applied perpendicular to the graphene. The Dirac 




FIG. 1: The magnetic field applied perpendicular to the 
graphene is zero in the center region of width 2W , and at the 
two sides of the strip they are in opposite directions with the 
same magnitude B. The total width of the strip is L > W. 
The classical trajectory of a typical snake state is also drawn 
schematically. 

Hamiltonian of graphene in low energy approximation 
reads H± — v (a x ir x ± <j y ir y ), where the index + (— ) 
corresponds to the valley K (K 7 ), the Fermi velocity 
is v — y/3/2aj/H given by the coupling constant 7 be- 
tween the nearest neighbours in graphene (a = 0.246 nm 
is the lattice constant in the honeycomb lattice), while 

= (tt x , n y ) = p — eA is given by the canonical momen- 
tum p and the vector potential A related to the magnetic 
field as B = rotA. The Pauli matrices a x and a y act in 
the pseudo-spin space. The Hamiltonian is valley de- 
generate for arbitrary magnetic field, ie, <j x H±<j x = H^, 
therefore it is enough to consider only one valley. In what 
follows we take valley K. 

The energy spectrum of our system can be obtained 
from the Schrddinger equation H + ty(x,y) = E^>(x,y). 
In our analytical calculations we assumed that L I, 
where I = y/H/\eB\ is the magnetic lenght. To construct 
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the wave function in each region we utilize the symme- 
tries of the system. The system is translation invariant in 
the y direction and [H, p y ] — in Landau gauge if the vec- 
tor potential has a form A = (0, A y (x), 0) T . Therefore 
the wave function can be separated: #(x, y) = <&{x)e %ky , 
where k is the wave number along the y direction. One 
can also show that [H, a y T x ] = 0, where T x is the reflec- 
tion operator, ie, it acts on an arbitrary function f(x) 
as T x f{x) — f{—x). This implies that the wave func- 
tions <&(x) can be classified as even wave functions satis- 
fying a y T x <&( e \x) — <f>( e '(x) and odd functions, when 
a y T x ^°\x) = -$ (o) (x). Moreover, it is true that 
[cr y T x ,p y ] = 0. From these commutation relations the 
even (odd) wave function ansatz in the three different 
regions can be constructed. In the central region, ie, for 
|x| < W and at energy E it is given by 
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where tany = k/K, K — \J e 2 — k 2 is the transverse 
wave number and e = E/(hv). In the left hand side, ie, 
for x < —W the wave function in Landau gauge reads 
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where £ = -y/2 (x + W + kl 2 ) /I, a± = (±1 - e 2 / 2 ) /2, 
7] = —iy/2/(el), and U(a,x) is a parabolic cylinder 
function (from the two parabolic cylinder functions 
we take that which tends to zero for x — > — oo). The 
wave function ansatz in the right region, ie, for x > W 



can be obtained from <J>^' o) (x) as <f>${x) = cr y T x ^'(x) 
and $^(x) = -a y T x ^ ( f\x). When k > e, the above 
given transverse wave number K has to be replaced by 
K = —iy/k 2 — e 2 . The boundary conditions require the 
continuity of the wave function at x = ±W resulting 
in a homogeneous system of equations for the ampli- 
tudes ci and C2- Hence, nontrivial solutions can be 
obtained from the secular equation resulting in energy 
bands E n (k) labelled by n = 0, ±1, ... for a given k. Ow- 
ing to the chiral symmetry (a z H±a z = —H±) we have 
E- n (k) =-E n {k)\llk. 

Figure [2] shows the comparison of energy bands E n (k) 
around the K point obtained from the Dirac equation 
and from tight binding (TB) model. Note that our TB 
calculation is similar to that made by Wakabayashi et al. 
except that they applied uniform magnetic field [17j. In 
the figure the surface states calculated from TB model 
are not shown. This will be discussed below. The ener- 
gies are scaled in units of Huj c , where tkj c = \f2hvjl — 
y/3/2'ja/l. One can see that the agreement between the 
two calculations is excellent. For large enough positive 
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FIG. 2: The energy bands E„(k) (in units of huu c ) as a func- 
tion of k (in units of 1/W) around the K point for magnetic 
field corresponding to W/l = 2.2. The solid lines (open cir- 
cles) are the results obtained from the Dirac equation (TB 
model, [l6]). Only the conduction band (E n (k) > 0) with 
n — 0,1, ... ,8 states is shown (n = corresponds to the low- 
est conduction subband, and the other bands are in increasing 
order in energy). The even (odd) wave function states corre- 
spond to even (odd) quantum number n. States Ai, Bi and 
Ci at energy 0.688 hu> c (dotted line) are snake states (see the 
text). 



k each state evolves into dispersionless, twofold degen- 
erate Landau levels having the same energies as in uni- 



, (lc) form magnetic field, ie, E^(k) = sgn(n) ftw c 



with 

n = 0, ±1, . . . , where sgn(-) is the sign function (see, eg, 



Ref. OH)- However, for negative wave number k the 
energy bands are dispersive. Examples for such states 
are A\ , B\ and C\ in Fig. [2l Each of these states carries 
current in the y or — y direction depending on their group 
velocity. The corresponding wave functions are localized 
in the central, zero magnetic field region as shown below 
(see Fig. [3]). These are the snake states in our system. 

It is also clear from Fig. [2] that the number of left and 
right moving states at a given energy are not the same 
which seems paradoxical at first sight. The ground state 
of the system appears to be unstable thermodynamically 
beacuse there is a net current flow in the —y direction 
even in equilibrium. To understand this paradox one 
needs to consider the surface states localized at the edges 
of the system (x = ±L/2). Figure [3] shows the same en- 
ergy bands E n (k) as in Fig. [2] obtained from TB calcula- 
tions for all allowed k. Two extra subbands appear and 
the resulting states at a given energy are denoted by D± 
and £>2 in Fig. [3] (although they are hardly distinguish- 
able because of the parameters we used, see the caption of 
Fig. [2]). We shall show that these non-degenerate states 
are surface states caused by the finiteness of the sample 
in the x direction and carry current in the y direction at 
the two edges, respectively. Moreover, in Fig. [3J states 
Ai,2, -Bi,2, and Ci,2 are snake states and the first two 
states carry current in the —y, while C±,2 in the y direc- 
tion. For finite L the surface states at the edges of the 
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FIG. 3: The dispersion relation E„(k) for all allowed k (in 
units of 1/a) including a small portion of the valence bands. 
The parameters are the same as in Fig. [2] for TB calculation. 
At the same energy as in Fig. [5] indicated by dashed line there 
are eight states: Ai, Bi, d and Di with i = 1, 2 (see the text). 



sample will compensate the net current carried by the 
snake states. It is easy to see that including the surface 
states the number of left and right moving states are the 
same for all energies ensuring a stable equilibrium state 
of the system. We expect the same result using the Dirac 
Hamiltonian for finite L (see works by Brey and Fertig, 
and Abanin et al. in Ref. 0]). 

To get better insight into the nature of the snake and 
surface states we calculated the current density distribu- 
tion of these states shown in Fig. 0J As can be seen from 
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FIG. 4: Current density distributions j y (x) (in arbitrary 
units) as a function of x (in units of lattice site in x rang- 
ing from 1 to N) for states A\, B\, C%, Di and D2 shown in 
Fig. E] 

the figure the snake states A\ and B\ carry current in 
the — y direction, while state C\ in the opposite direc- 
tion. Moreover, all of these states are localized in the 
center of the sample. States D\ (-D2) are surface states 
with current flowing in the y direction close to the left 



(right) edges of the sample. States Ai, Bi and Ci with 
i = 2 behave the same way as that with i = 1 due to 
valley degeneracy. 

Varying the Fermi energy Ep the character of these 
states would change and consequently their current dis- 
tribution too. For Fermi energy lying between the first 
and second Landau levels, ie, for E^k) < Ef < E^{k) 
the net current contributing from snake states B\ and 
Ci localized at the center of the sample is zero, therefore 
they are locally compensated states. While the snake 
state Ai is also localized at the center, it is locally un- 
compensated. Only the current from surface state Di 
will compensate the current from the snake state A\ to 
ensure the stability of the ground state of the system. 
However, locally this snake state is uncompensated since 
the overlap between states A\ and D\ is negligible. When 
the Fermi energy crosses the Landau level n = 1 then not 
just state Ai but B\ also becomes uncompensated, while 
state C\ evolves into a surface state. Thus, we find that 
in graphene for all Fermi energies there is always at least 
one locally uncompensated snake state localized at the 
center of the sample and it is compensated only with 
surface states far from the center part of the system. 

It is instructive to compare these results with that ob- 
tained for 2DEG with the same magnetic field profile as 
in Fig [T] As can be seen from Fig. [5] the Landau levels 
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FIG. 5: The energy bands E n (k) of 2DEG obtained from TB 
model for the same magnetic field profile as in Fig. [T] and with 
magnetic strength as in Fig. [2] The energies are in units of 
hQ. c , where Q c = eB/m with effective mass m of electrons, 
and k is in units of 1/a, where a is the lattice constant in the 
square lattice). Here W = L/10 and L = 200a. 

are at energies E n (k) = HD, c (n + l/2), where n = 0, 1, . . . 
as expected in case of uniform magnetic field. Using a 
similar analysis as before for graphene one can find snake 
and surface states 19]. The dispersion relation implies 
that in 2DEG no locally uncompensated current appears 
below the first Landau level (n = 0). This is a striking 
difference between the two systems. 

In summary, we studied the dynamics of electrons in 
graphene when the applied magnetic field is non-uniform. 
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For a simple, step-like magnetic field dependence we show 
that snake states appear similarly to the case of 2DEG. 
We calculated the energy bands in case of an infinite sys- 
tem using the Dirac Hamiltonian which agree very well 
with that obtained from our TB calculations. Moreover, 
we show that the surface states of a finite width sample 
ensure the stability of the ground state of the system. 
We find that in contrast to 2DEG the snake states in 
graphene are locally uncompensated for all Fermi en- 
ergies. The different behavior of the snake states in 
graphene compared to the 2DEG is a clear manifestation 
of the massless Dirac fermion like excitations in graphene. 
The current carrying snake state at low enough Fermi 
energy (between the Landau level n = and n = 1 or 
n = —1) is expected to be as resilient against scattering 
on impurities as the surface states in quantum Hall effect 
of graphene. The peculiar behavior of the snake states 
in graphene with non-uniform magnetic fields could be 
utilized in future experimental and theoretical works. 
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